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A two-flavor color superconducting (2SC) Nambu-Jona-Lasinio (NJL) model is introduced 
at finite temperature T , chemical potential /i and in the presence of a constant magnetic 
field eB. The effect of (T, /i, eB) on the formation of chiral and color symmetry breaking 
condensates is studied. The complete phase portrait of the model in T — — eB, and 
T — eB phase spaces for various fixed eB, T , and is explored. A threshold magnetic field 
eBt ~ 0.5 GeV^ is found above which the dynamics of the system is solely dominated by the 
lowest Landau level (LLL) and the effects of T and are partly compensated by eB. 
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I. INTRODUCTION 

Recently, the study of the properties of quark matter in the presence of strong uniform magnetic 
fields has attracted much attention. Possible effects caused by strong magnetic fields include magnetic 
catalysis , modification of the nature of electroweak phase transition [3], spontaneous creation 
of axial currents [J], formation of ttq domain walls and chiral density waves ^, chiral magnetic 
effect 0,0], and last but not least the influence on possible color-superconducting phases [s-lG]. In 



this paper, we will focus on the magnetic catalysis and its possible effects on the phase diagram 
of a magnetized two-flavor color superconducting NJL model at finite temperature and chemical 
potential.^ In particular, the dependence of the included meson and diquark masses on thermodynamic 
parameters, and possible interplay between these parameters on the formation of meson and diquark 
condensates and on the nature of phase transition will be scrutinized. 

At zero temperature, it is known that strong magnetic fields enhance the production of chiral and 
diquark condensates, albeit through different mechanisms, as it is described in Whereas magnetic 
catalysis of chiral symmetry breaking [l| is mainly responsible for dynamical mass generation and en- 
hances the production of chiral condensates by increasing the particle- antiparticle interaction strength, 
a certain modification in the density of states of charged quarks near the Fermi surface, depending on 
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the external magnetic field, reinforces the pairing of charged quarks and is made responsible for the 
enhancement of diquark production by a penetrating strong magnetic field In other words, in 



contrast to the effect of magnetic catalysis of chiral symmetry breaking, which is essentially based on 
a dimensional reduction of the dynamics of fermions from Z) = 3+ ltoD = l + l dimensions to the 
regime of LLL dominance [l| , the pairing mechanism by color superconductivity, involving the charged 
quarks near the Fermi surface, is already D = 1 + 1 dimensional, and therefore the external magnetic 
field does not lead to any further dimensional reduction It is the goal of the present paper to 



explore possible effects of finite T and fi on the above mechanisms of chiral and diquark production 
in the presence of strong magnetic fields. We will show that in our setup a certain threshold magnetic 
field exists, above which the dynamics of the system is solely dominated by the lowest Landau level 
and the effects of T and /x are partly compensated by very large eB. The largest observed magnetic 
field in nature is ~ 10^^ — 10^'^ Gauf5 in pulsars and up to ~ 10^^ — 10^^ GauB on the surface of some 
magnetars, where the inner field is estimated to be of order ~ 10^^ — 10^^ Gauf^ [22]. In the early 



23]. Su- 



universe, magnetic fields of order ~ 10^'' Gaui5 may be produced at the beginning of infiation 
perconductiye cosmic strings, if they exist, may have magnetic fields up to ~ 10^^ — 10^^ GauB in their 
vicinities 



241 ]. In 25|], it is shown that "the maximum value of magnetic field that delimits the range 
of values admitted without revising QED" is of order ~ 10^^ GauB. There are also evidences for strong 
magnetic field creation in non-central heavy ion experiments |26l ]. The early estimate of the magnetic 



field for the RHIC energy was made in |27i], where it was shown that the magnitude of the magnetic 
field for an earlier stage of noncentral Au-Au collision at energy ^snn = 200 GeV and the impact 
parameter ~ 4 fm is about eB ~ 1.3 ~ 0.025 GeVf, that corresponds to ~ 4.3 x 10^^ Gaufi.^ 
Using a microscopic transport model, the authors in |28|] estimate the lowest bound of the maximal 
magnetic field strength at the LHC energy ^/snn = 4.5 TeV with the same impact parameter, to be 
eB ~ 15 ~ 0.3 GeV^, which is equivalent to -B ~ 5 x 10^^ GauB. Numerically, it is not a priori 
clea. „he.he, th.. amounts of .agnet.c fields „e .a.ge eno.gh to ,us..ty .he LLL app_on, 
as it is done in [1] to demonstrate the magnetic catalysis of chiral symmetry breaking and in [11[ to 
demonstrate the effect of strong magnetic field on color superconductivity. 

We have tried to answer to this question, among others, in [l6], where a color neutral and dense 
2SC-NJL model has been introduced in a (rotated) constant magnetic field at zero temperature. To 
do this, first the dependence of chiral and color symmetry breaking condensates on the chemical 
potential /x and the rotated magnetic field eB is determined analytically in a lowest Landau level 
(LLL) approximation. Then, the meson and diquark masses are computed numerically for arbitrary 



Here, = 140 MeV. Moreover, eB = 1 GeV^ corresponds to B = 1.69 x 10^" Gaufi. 
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magnetic fields. Comparing these analytical and numerical data, we have found a certain threshold 
magnetic field eBt — 0.45 — 0.5 GeV^, corresponding to ~ 8 — 8.5 x 10^^ Gaufi, above which the 
system turns out to be dominated by LLL. Below eBt, the chiral and diquark mass gaps oscillate with 
the external magnetic field. These oscillations are the results of the well-known van Alfven-de Haas 
effect 



m 



21 



291. that occurs whenever Landau levels pass the quark Fermi surface. They are also observed 
30| , where the dependence of chiral symmetry breaking mass gaps on constant magnetic fields is 



explored. Similar oscillations are also perceived in 



14l ] , where diquark mass gap and the magnetization 



corresponding to the superconducting magnetized color-flavor locked (CFL) phase are determined as 
a function of external magnetic fields. For eB > eBt, the system enters a "linear regime", where 
the mass gaps and the magnetization depend linearly on the external magnetic field. In this regime, 
where the system is believed to be solely affected by the dynamics of the fermions in the LLL, the 
numerical data coincide with the ana,lytical results (see [l6| for more detail). What concerns the phase 



diagram of the model introduced in [16|], it is shown that a first order phase transition occurs between 
the chiral symmetry breaking (xSB) and the color symmetry breaking (CSC) phases at T = and 
H ~ 350 — 450 MeV. This transition is then followed by a second order phase transition between the 
CSC phase into the normal quark matter.^ 

In the present paper, the same magnetized 2SC model will be considered at finite temperature. We 
are in particular interested on the additional effects of finite temperature, and, will focus on possible 
interplay between T, /x and eB on the formation of chiral and color symmetry breaking condensates and 
on the nature of phase transitions. Our results may be relevant for the physics of heavy ion collisions, 
where, recently the question of accessibility of color superconducting quark matter phases is pointed 



out 3^]. The authors in [3^ use a Polyakov-NJL (PNJL) model at finite temperature and density 
and present the corresponding QCD phase diagram including, among others, mixed phase regions of 
first order transition of 2SC-CFL quark matter and second order 2SC-normal phase transition. From 
this phase diagram, they conclude that color superconductor phase is already accessible at the present 
nuclotron-M energies 4 < < 8 AGeV, and, that possible transition from 2SC to normal quark 
matter becomes attainable in the planned FAIR-CBM and NICA-MPD experiments at 2 < ii^ < 40 
AGeV. In this paper, the additional effect of constant magnetic field, that is believed to be created 
in non-central heavy ion collisions will be explored from purely theoretical point of view. A complete 
phenomenological answer to the question of accessibility of color superconducting phases in heavy ion 



^ As it is known from |31l |. in the regime of low temperature and large chemical potential, the 2SC phase goes over into 
the three-flavor CFL phase. In a two- flavor model, however, where no CFL phase can be built, only a simple transition 
from the 2SC to the normal phase is assumed to exist. Same assumption is also made in As we will note in Sec. 
Ill, our numerical results for low temperature and large chemical potential will be only of theoretical nature. 
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collisions is out of the scope of this paper, and will be postponed to future publications. 

The organization of the paper is as follows: In Sec. II, we introduce our model for the hot mag- 
netized two-flavor superconducting quark matter and determine the corresponding thermodynamic 
potential. In Sec. Ill, our numerical results on the dependence of the meson and diquark gaps on 
chemical potential, magnetic field, and temperature will be presented in Sees. III.A.l - III. A. 3. In 
Sec. III.B, the complete T — fi, T — eB and fi — eB phase portraits of the system at various fixed eB, fi 
and T will be illustrated. The above mentioned threshold magnetic field, eBt, will be determined by 
comparing the analytical and the numerical results, corresponding to the second order critical lines of 
the transition from the xSB and CSC phases to the normal phase, in a T — plane for eB = 0.5, 0.7 
GeV^ (Sec. III.B. 1) and in aT-eB plane for /x = MeV (Sec. III.B. 2). The details of the analytical 
computations, that lead to second order critical surfaces of these transitions in a (T, ^, eB) phase 
space, will be presented in App. A. Section IV is devoted to concluding remarks. 

II. MAGNETIZED 2SC QUARK MATTER AT FINITE T AND n 

In Sec. II. A, we briefly review our results from [l6^ and introduce a two-flavor NJL model including the 
meson and diquark condensates at finite temperature and density and in the presence of a constant and 
uniform (rotated) magnetic field. In Sec. II. B, the corresponding one-loop thermodynamic potential 
will be determined in a mean field approximation. 

A. The Model 

The Lagrangian density of a two-flavor gauged NJL model is given by 

£f = ?(x)[i7'^(a^ - ieQA^ - igT^G^) + ^l7'>]^^^{x) 

+Gsmxmx)f + (V^(x)i75rV'(x))2] + GD[{ii^^{x)efeh5i^{x)){ii^{x)efeh5i;^{x))]. (II.l) 

Here, = Cip and tp = ip'^C are charge-conjugate spinors, and C = ^7^7*^ is charge-conjugation 
matrix, r = {ti,T2,t^) are Pauli matrices. Moreover, {£f)ij and (ec)"'^ = (^c)"*^ are antisymmetric 
matrices in color and flavor spaces, respectively. For a theory with two quark flavors and three color 
degrees of freedom, the flavor indices i,j = (1,2) = {u,d), and the color indices a,b = (1,2,3) = 
{r,g,b). The quarks are taken to be massless = = 0. The quark chemical potential which 
is responsible for nonzero baryonic density of quark matter is denoted by /x. Here, = where 
= diag(l,l,— 2) the 8*^ Gell-Mann A-matrix. The scalar and diquark couplings are denoted 
by Gs and Gd, respectively. The charge matrix Q = Qj Ic, where Qf = diag (2/3, — 1/3) is the 
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fermionic charge matrix coupled to U{1) gauge field A^. The same setup with an additional color 
chemical potential /ig, imposing the color neutrality of the theory, is also used in [IS] to study the 



effect of magnetic field on quark matter under extreme conditions. Following the same steps as in jl6l ] 
to determine the Lagrangian density containing the chiral and diquark condensates in an appropriate 
Nambu-Gorkov form, we define first the mesonic fields 

a = -2Gs(iptl^), and n = -2Gs{^i-r'^T^), (II.2) 

as well as the diquarks fields 

A = -2GD{ii^^ef€h5i^), and A* = -IGoiii^Efeh^i^^). (II.3) 

Combining then the gauge fields and G^, using the "rotated" charge operator Q = Qf®lc — 
I —2= J J the rotated massless Uemi^) field, A^ = A^ cos sin 9, as well as the massive in- medium 

S*'' gluon field, G^ = A^sinO + G^cosO can be derived (see [161 ] for more details). Replacing A^ with 
an external gauge field A^^* = {0,0, Bx,0) in the Landau gauge, a constant rotated background U{1) 
magnetic field directed in the third direction B = Be-^ is induced. Neglecting then the massive gauge 
boson G^, the total modified bosonized Lagrangian density, C = Ck + Cf, in the presence of a uniform 
magnetic field arises. It consists of a kinetic term 

and an interaction term 

(11.5) 

Assuming that the vacuum of the system is characterized by (cr) 7^ and (vf) = 0, we have neglected 
the 7? mesons. Moreover, using the definition of the rotated charge operator Q in a 6-dimensional 
flavor-color representation {ur,Ug,Ub,dr,dg,db), the rotated q charges of different quarks, in units of 
e = e cos 6, are given by 



quarks 


Ur Ug Uf, dr dg dh 


q 





To bring the above Lagrangian density £j in a more appropriate Nambu-Gorkov form, we introduce 
at this stage the rotated charge projectors 

Oo = diag(0, 0, 0, 0, 0, 1), Oi = diag(0, 0, 1, 0, 0, 0), 

n^i = diag(l, 1,0, 0, 0, 0), 0_ 1 = diag(0, 0, 0, 1, 1,0), (II.6) 



that satisfy QQg = qytq. The Nambu-Gorkov bispinors are then defined by 



where ipq{x) = Qqip{x). In terms of ^'g and ^q, the Lagrangian density Cf from pi.5P in the Nambu- 
Gorkov form reads (for more details see [3]) 

9e{o,i,±i} 

For q G {0, 1}, Sq is given by 
and for q G { — ^, +|}; reads 



[G7,J-i 



Here, 

[Gf-)]-i = 7''(^9^ + - ^ ± A^'^/^o), (11.10) 

and K^^"^ = ^At2"'A2^7^^ as well as k' = ■Jqk^Jo = iA*r2A27^. In the next section, the Lagrangian 
density £ with £fc in pi.4p and £/ in pi.7p - pi.l0p will be used to determine the thermodynamic 
potential of this model in the mean field approximation. 



B. Thermodynamic potential 

The quantum effective action of the theory, T^g, is defined by integrating out the fermionic degrees of 
freedom using the path integral 



gir,ft[<x,A,A*] ^ J p^p^exp (^i j d\ £^ . (11.11) 



At one loop level, it consists of two parts: the tree-level and the one-loop effective action, and 
^[l^K In the mean field approximation, where the order parameter a = ((t(x)), A = (A(a;)) and 
A* = (A*(x)) are constant, the tree-level part of F^ff is given by 

rS)[,.A.A-;B,.-v(|-.l|^.f). (ILO, 
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Here, V is a four-dimensional space-time volume and |Ap = AA*. The one- loop effective action is 
given by 

Til\a,A,A*;B] = -l^TV^^,^,, HSg'^], (11.13) 

where Sq is defined in pi.8P - (|II.10p . In (jll.lSP the trace operation (Tr) includes apart from a two- 
dimensional trace in the Nambu-Gorkov space (NG), a trace over color (c), flavor (/), and spinor 
(s) degrees of freedom as well as a trace over a four-dimensional space-time coordinate (x). After 
performing the trace operation over the NG,c, /, and s using the method described in we arrive 
at 

fliHp-)= E f«/^(p-), (11.14) 
that includes the contribution of the blue (6), red (r), and green (g) quarks 

ge{o,i} 

E flfi^^'(p) = -2i E lndet,.[(4+^)'-pg)(4-^)'-pg)]. (11.15) 



Here, p is a modified four- momentum defined by [s] 



(n.l6) 



P'^^o = (Po,0, ill Y^2|geS|n,p3), for q = l,±i, 

P^=o = (Po,p), for g = 0, 

with p = {pi,P2,P3)- 111 (III.15p . Eg are given by the dispersion relations corresponding to the neutral 
and charged particles [ij] 

Eg = ^2\qeB\n + pi + a^ for q = l,±i, ^^^^^^ 

Eo = Vp2 + a2, for = 0. 

Moreover, we have E^^^ = ^ {Eg it /i) + | A^. Performing the remaining determinant in the coordi- 
nate space, a space-time volume factor V arises. Combining then the resulting expression with the 
tree- level part of the effective action from 01. 12^ . the effective action of the theory can be given in 
terms of the effective thermodynamic (mean field) potential ^l^g as T^g = — V$7off. Introducing now 
discrete Matsubara frequencies by replacing pQ with pQ = iuj£, where uj£ = ^{2£ + 1) and /3 = T^^, the 
one-loop effective potential = —ypT^g^ is first given by'* 

p e=-oo K=±i ^ gefoA} J 



(n.l8) 



^ In imaginary time formulation, the four-dimensional space-time volume V is replaced by V ^ VP, where V is the 
three-dimensional space volume and /? — 1/T the compactification radius of imaginary time coordinate. 
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Converting then the logarithms into proper-time integrals over the dimensionful variable s, and using 
the Poisson resummation formula 



oo 



y =^[i + 2ycos(x/3^)e-^], (11.19) 



£=-oo ^ e=i 

to separate the resulting expression into a temperature dependent and a temperature independent 
part, we arrive at the one-loop effective potential of our model 



oo 2 2 

x[l + 2j;(-l/e-^], (11.20) 

where -Eg are defined in ()II.17p . Replacing at this stage the discrete sum over momenta with continuous 
integrations over momenta by making use of 



for neutral, and 



+°o „+oo 
p 71=0 

for charged particles,^ and eventually adding the tree level part of the effective potential to the resulting 



1 r+oo J 

-^f{P^^o)^\qeB\Y.anJ ^ ^ f{n,P3), (11.22) 



expression, the full mean field effective potential at one-loop level is given by^ 

"^ ' '^^^"4^5 + 40^,+ 2 +2V^^^jy (2vr)3X 

+eB±a^ r^(e-(^-+-)^ +2e-W^l%,)|[l + 2f;(-l)^ e'^^]. (11.23) 

Note that in ()II.22p as well as pi.23p . n denotes the discrete Landau levels and an = 2 — Sno is 



14l |. Moreover, 



introduced to consider the fact that Landau levels with n > are doubly degenerate 
we have used E+1/2 = ^-1/2- I^i the next section, we will use pL23p to determine numerically the 
chiral and diquark gaps and to present the complete phase structure of the magnetized two-flavor 
superconducting NJL model at finite T and /i. 



^ Only charged particles interact with the external magnetic field. 
In [16| . the same effective potential (|II.23I) was determined using a different method (see (3.19)-(3.22) in 
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III. NUMERICAL RESULTS 



In the previous section, we have determined the effective potential (|II.23p of two-flavor NJL model 
including meson and diquark condensates at finite temperature, chemical potential and in the presence 
of constant magnetic fields in the mean field approximation at one-loop level. It is the purpose of 
this paper to have a complete understanding on the effect of these external parameters on the quark 
matter in the 2SC phase. This will complete our analysis in where only the effect of fi,eB was 
considered at T = 0. We start this section with presenting the numerical results on the fi, T and eB 
dependence of the chiral and diquark condensates. We then continue with exploring the T — /j,, T — eB, 
and fi — eB phase diagrams for fixed values of eB, fi and T, respectively. Before presenting our results, 
we will fix, in the subsequent paragraphs, our notations and describe our numerical method. 

To determine the chiral and diquark gaps, the thermodynamic potential O^ff from pi.23p is to be 
minimized. To solve the corresponding gap equations 



dn,a{^,A;T,fi,eB) 



9J^ia,A;T,,,eB) 



rr A dA 



= 0, (III.l) 

o-bAb 



da 

numerically, we have to fix the free parameters of the model. Our specific choice of parameters is 



16|, 



33| 



A = 0.6533 GeV, Gs = 5.0163 GeV-^, and Go = ^Gs, (III.2) 

where A is the momentum cutoff and Gs as well as Gd are the chiral and diquark couplings. Using 
this special set of parameters, we have shown in [l^, that at T = no mixed phase characterized, 
with (as 7^ 0,Ab ^ 0), will appear. The same feature persists at finite T. Moreover, for vanishing 
magnetic field and at zero temperature, the parameters in ()III.2p yield the meson mass ctq — 323.8 
MeV at /X = 250 MeV, and the diquark mass Aq ~ 126 MeV at n = 460 MeV jl6].'^ To perform the 



momentum integration over p and in (|II.23P numerically, we have introduced, as in 16|], smooth 
cutoff functions (form factors) 

/a = ^A.^ = / , , X . (ni.3) 

l + exp(M^) l + exp(VS±ap^^ 

that correspond to neutral and charged particles, respectively.^ In pil.3p . j4 is a free parameter and is 
chosen to be ^ = 0.05A. Similar smooth cutoff function (form factor) is also used in jl^. Here, as in 



Although our free parameters A,Gd, and Gs coincide with the parameters used in [33 |. the numerical value of ao 
is different from what is reported in (3^. The reason for this difference is apparently in the choice of the cutoff 
function. Whereas in [sH a sharp momentum cutoff is used, we have used smooth cutoff function (|III.3P to perform 
the momentum integrations numerically. 
* In pi.23|l , the integrals proportional to eB and including a summation over Landau levels n arise from charged quarks 
with charges q = ±i, +1. 
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[13|, the free parameter A determines the sharpness of the cutoff scheme. At this stage, let us notice 
that the solutions of pil.ip are in general "local" minima of the theory. Keeping (a, A) 7^ (0, 0) and 
looking for global minima of the system described by complete f^cft^c, A; T, fi, eB) from pi.23p . it turns 
out that in the regime /x e [0, 800] MeV, T £ [0, 250] MeV and eB G [0, 0.8] GeV^, the system exhibits 
two "global " minima. They are given by (as 7^ 0,Ab = 0) and (as = 0,Ab ^ 0). We will denote 
the regime characterized by these two global minima by the xSB and the CSC phases, respectively. 
According to the above descriptions, in order to determine the chiral and diquark condensates, we will 
use instead of the gap equations (llll.lj) . 

dn,s{a,AB = 0;T,fi,eB) 



da 



(III.4) 



in the xSB phase and 



dn,^{aB = 0,A;T,fi,eB) 



dA 



0, 



(III.5) 



in the CSC phase to simplify our computations. As it turns out, apart from the xSB and the CSC 
phases, there is also a normal phase characterized by {ctb = 0,Ab = 0). Following the method 



presented in 



20|, we will determine the critical lines of first order transition between the xSB phase 



and the normal phase using 
Ab = 0) 



da 



and 



aft (ctb, Ab = 0) = aft (ctb = 0, Ab = 0). 



(III.6) 



Similarly, the first order phase transition between xSB and the CSC phases is determined by solving 



aaff(c7, Ab = 0) 



da 



dn^sicTB = 0, A) 



OA 



0, 



(III.7) 



and 



aff(cTB = 0, Ab) = aff(c7B, Ab = 0), 



(III.^ 



simultaneously [20]. The second order critical lines between the xSB as well as CSC phase and the 



normal phase will be also determined by the method described in 



20( (see Eq. (2.35) in [20]): To 



determine the second line between the xSB and the normal phase, we solve 

aaff(c7,A = o) 

hm —5 = 0. 

(j2^o oa^ 

The second line between the CSC and the normal phase is then determined by solving 

r gaft(cT = 0,A) 

hm ——^ = 0. 

A2^o dA^ 



(III.9) 



(III. 10) 
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The same method is also used in [2 ill to find the second order phase transition between the xSB and 
the normal phase (see page 9 in 2l[). To make sure that after the second order phase transition the 
global minima of the effective potential are shifted to cr = in pil.op and to A = in (jlll.iop . and in 
order to avoid instabilities, an analysis similar to [ol is also performed. Same method is also used in 



[161 ] to determine second order critical line between the CSC and the normal phase. 



A. The (/i, T, eB) dependence of chiral and diquark condensates 



1. The pL- dependence of ub and As 



The ^-dependence of both gaps at different temperatures, T = 0, 20, 70, 150 MeV, and magnetic 
fields, eB = 0,0.3,0.5 GeV^, are demonstrated in Figs, [m panels (a)-(c), respectively. The green 
dashed and solid lines denote the aB mass gaps. The red solid lines determine the normal phase with 
aB = As = 0. The diquark gaps A^ are demonstrated with blue solid lines. Dashed (solid) lines 
denote the first (second) order phase transitions. 



a) T = MeV, eB = GeV^ 
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FIG. 1: The /^-dependence of ctb and As is demonstrated for T = MeV and eB = 0, 0.3, 0.5 GeV^ in (a), (b) 
and (c), respectively. The transitions from xSB to the CSC phase is of first order, the transitions from the CSC 
to the normal phase is second order. 

The plots in Figs. [Ua)-[T][^c) show that at T = MeV, the magnetic field enhances the formation of 
chiral condensate ctb- The value of aB is constant in /i, with fi < fic and fic — 320 — 350 MeV, and 
increases with eB. On the other hand, for small value of eB < 0.5 GeV^, Ab increases with in 
the regime fic < fJ- < 600 MeV [Figs. (H^a) andlTJb)]. Similar observation is also made in [l^, where 
the same model as in the present paper is studied, and additionally the color chemical potential /ig is 
assumed to be nonzero. Comparing the diagrams of Fig. [T] with the corresponding diagrams in Fig. 
7 of jiol, it turns out that fis has no significant effect on the /^-dependence of fj^ and A^ at T = 
MeV.^ For eB = GeV^, the /i-dependence of diquark mass gap in Fig. [U^a) can be compared with 



In [35|], in eB = case, the color chemical potential /ig is shown to be small and its effect is therefore neglected. See 
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the analytical result 



^0 



C2(A^ - ^^)exp (^-^ - V j ' 

from In (IIII.llI) . C2 = 46"^ ~ 0.2 and gd = ^^5^. Plugging the numerical value of A and 
from piL2p in pil.lip and plotting the resulting expression in a Aq vs. fi diagram, the result is in 
close agreement with the /i-dependence of diquark mass gap in Fig. H^a) in the regime 320 < fi < 600 
MeV. In the same regime of /U, the //-dependence of Aq from ()III.11|) agrees also with the well-known 
34| . In contrast, as it can be seen in Fig. H^c), for eB = 0.5 GeV^, A^ decreases with 



result from 



fi G [350, 600[ MeV. This behavior is expected from the analytical expression 



A2 



that is also computed in [3] using an appropriate LLL approximation.^'^ In piI.12|l . Ab = VaB and 
. Plugging the numerical values of A and Gd from (jTTL2]) in ()III.12p . and plotting the 



4(Al,-M=)exp|-£l 



(III.12) 
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resulting expression in a A^ vs. /j, diagram for eB = 0.5 GeV^, it turns out that A^ decreases with 
^ in the regime 350 < n < 600 MeV, as it is shown in Fig. [D^c). Similar expression as ()III.12|) is 



also derived in 



111 ] for the diquark mass gap in a three-flavor CFL model in the presence of strong 



magnetic field, using the LLL approximation. 

In Figs. [2ja)-l2jc) the /^-dependence of as and A^ is plotted at T = 20 MeV and for eB = 0, 0.3, 0.5 
GeV^. Comparing with the results from Figs. [Ha)-[l]J^c), it turns out that increasing temperature up 
to T = 20 MeV, has no significant effects on the results of T = MeV. This is in contrast with the 
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b) T = 20 MeV, eB = 0.3 GeV^ 
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FIG. 2: The /^-dependence of as and Ab is demonstrated for T = 20 MeV and eB = 0,0.3,0.5 GeV^ in (a), 
(b) and (c), respectively. The transitions from xSB to the CSC phase is of first order, the transitions from the 
CSC to the normal phase is second order. 

situation at T = 70 MeV. In Figs. [3l^a)-l3jc), the /^-dependence of aB and A^ are plotted at T = 70 



also [3g|, where two-flavor magnetized color superconducting quark matter with /ig = is studied. 
Later we will see that eB ~ 0.5 GeV^ is strong enough to justify LLL approximation. 
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MeV and for eB = 0,0.3,0.5 GeV^. As it turns out, in the regime fj, < fic, the chiral condensate as 
increases with the magnetic field. The diquark condensate appears in the regime fi G [480,600[ MeV 
for eB = 0,0.3 GeV^ [Figs. EJa) andEJb)]. For eB = 0.5 GeV^ in Fig. EJc), however, no diquark 
condensate appears in the relevant regime n < 600 MeV. Moreover, comparing Figs. [2] and El we 
notice that at T = 70 MeV, in Fig. [3l in contrast to the situation at T = 20 MeV, in Fig. [2l the 
first order transition from the xSB to the normal phase does not occur over the CSC phase; in Fig. 
El for eB < 0.5 GeV^, there is a first order phase transition from the xSB to the normal phase, then 
a second order phase transition occurs from the normal to the CSC phase. For eB = 0.5 GeV^, the 
phase transition from the x^B to the normal phase is of first order but no CSC phase appears in the 
regime /x < 600 MeV. 



a) T = 70 MeV, eB = GeT 
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b) T = 70 MeV, eB = 0.3 GeV^ 
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c) T = 70 MeV, eB = 0.5 GeV^ 
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FIG. 3: The /i-dependence of as and As is demonstrated for T = 70 MeV and eB — 0, 0.3 GeV^ in (a) and (b), 
respectively. While the transitions from x^B to the normal phase (solid red line) is of first order, the transitions 
from the normal to the CSC phase is of second order, c) The /^-dependence of as is demonstrated for T = 70 
MeV and eB = 0.5 GeV^. The transition from the xSB to the normal phase is of first order and no CSC phase 
appears at /i > 330 MeV. 



At higher temperature, as it is demonstrated for T = 150 MeV in Fig. [H no diquark condensate 
appears at all in the relevant regime fi < 600 MeV [Figs. [l|^a)-[3|^c)]. Moreover, whereas for eB = 0, 0.3 
GeV^, the transitions from the xSB to the normal phase are of second order, for eB = 0.5 GeV^, the 
xSB-normal phase transition turns out to be of first order. The above observations suggest that while 
the magnetic field catalyzes the formation of the chiral condensate as, the diquark condensate is 
suppressed in the presence of a constant magnetic field. However, as it turns out, the latter feature 
depends on the strength of the magnetic field. This will be shown explicitly in the next section, where 
the effect of arbitrary magnetic fields on the formation of ctb and A^ for a wide range of eB G [0, 0.8] 
GeV^ will be explored. 
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a)T = 150 MeV, eB = OGeV^ 
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FIG. 4: The /^-dependence oi as is demonstrated for T = 150 MeV and eB = 0,0.3,0.5 GeV^ in (a)-(c). While 
for eB — 0, 0.3 GeV^ the transition from xSB to the normal phase is of second order, for eB — 0.5 GeV^ this 
transition is of first order. No CSC phase appears at T = 150 MeV. The solid red line denotes the normal phase 
with (ctb 0, As = 0). 

2. The eB -dependence of ctb and As 



The effect of external magnetic field on the formation of chiral condensates is studied intensively in 



the literature (see e.g. {l, Il7l4l9t|). As it is shown in [1], in the presence of strong magnetic field, 
the dynamics of the system is fully described by its dynamics in the LLL. However, there are, to the 
best of our knowledge, no evidences in the literature that fix quantitatively the strength of magnetic 



16l |. we have answered to this question 



field which is enough to justify a LLL approximation. In 
numerically. First, we have determined the eB dependence of the gaps in a LLL approximation, and 
then, compared the analytical results with the numerical ones including the effect of all Landau levels. 
It turned out that by increasing the magnetic field from eB = GeV^ to eB = 0.8 GeV^, first in 
the regime eB < 0.45 GeV'^, the mass gaps underly small van Alfven-de Haas oscillations. Above a 
certain threshold magnetic field eBt — 0.45 — 0.50 GeV^, the dependence of the condensates of eB was 
linear. In this linear regime, our numerical results were comparable with the analytical results arising 
from the solution of the corresponding gap equations in the LLL approximation. 

In Fig. O the eB dependence of meson and diquark mass gaps, as and A^, are plotted at T = 
and fi = 250 MeV [Fig. [5]^a)], and fi = 460 MeV [Fig. [^b)]. A comparison with similar plots from 
Fig. 2 of [l^ shows that color neutrality condition has no significant effect on the above mentioned 
threshold magnetic field eBt — 0.5 GeV^, and on the behavior of ub and below and above eBt, 
consisting of van Alfven-de Haas oscillations in eB < eBt and the linear rise in eB > eBt regimes. 
Similar strong van Alphen-de Haas oscillation of 2SC mass gap A^ in the regime eB G [0.4, 0.6] GeV^ 
in Fig. [5]^b) is also observed in jisl . Il4 | in the three-fiavor CFL phase in the presence of magnetic 



Sco Footnote 9. 
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a)fi = 250 MeV 



b)ii = 460 MeV 
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FIG. 5: (a) The dependence of the chiral gap as on eB for fixed n = 250 MeV, at T = MeV. (b) The 
dependence of diquark gap As on eB for fixed ^ — 460 MeV, at T = MeV. 



fields, albeit in another regime of eB/ fj?. 
a) 11 = 250 MeV 
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FIG. 6: (a) The dependence of the chiral gap ctb on eB for fixed /i — 250 MeV, at various temperatures, (b) 
The dependence of diquark gap As on eB for fixed /i = 460 MeV, at various temperatures. 

In Fig. Ella), the eB dependence of cr^ is demonstrated for fixed fi = 250 MeV and at various fixed 
temperatures T = 20, 60, 70 MeV; after small oscillations in the regime below the threshold magnetic 
field eBt, the system enters the linear regime, where only the contribution of the LLL affects the 
dynamics of the system. For eB < eBt, o'B decreases by increasing the temperature. In the linear 
regime eB > eBt, however, the effect of temperature is minimized. In Fig. [6lb), the dependence of the 
diquark gap is demonstrated for fixed ^ = 460 MeV and at various temperatures T = 20,60,70 
MeV. At T = 20 MeV, small oscillations occur in eB < eBt. In contrast, A^ monotonically increases 
with eB in the linear regime eB > eBt. The behavior of A^ at T = 20 MeV is not too much different 
than its behavior at T = MeV [see Fig. [5jb)]. This confirms our previous observation from the 
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comparison of Figs. [T]and[2J By increasing the temperature to T = 60 MeV, decreases. It vanishes 
in the range eB S [0.4, 0.6] GeV^. Comparing to the T = MeV plot from Fig. [5{b), it turns out that 
the value of at T = MeV in this regime of eB G [0.4, 0.6] GeV^ is in the same order of magnitude 
as 60 MeV. This is related with the well-known BCS ratio of critical temperature to zero-temperature 
gap in zero magnetic field case [33] [see the discussion below]. For eB > 0.6 GeV^, A^ increases again 
with eB. The situation demonstrated in Fig. [6]^b) at T = 60 MeV can be interpreted as a second 
order phase transition from the CSC to the normal phase at T = 60 MeV, fi = 460 MeV and eB ~ 0.4 
GeV^, followed up with a second order phase transition from the normal to the CSC phase at T = 60 
MeV, |U = 460 MeV and eB ~ 0.6 GeV^. We believe that this CSC-Normal-CSC phase transition are 
caused by strong van Alfven-de Haas oscillations in the regime eB G [0.4,0.6] GeV^. A comparison 
between the T = 60 MeV plot of Fig. EJ^b) and the fi = 460 MeV plot of Fig. [T^i) confirms these 
statements. Let us notice that van Alfven-de Haas oscillations, which are also observed in metals 



29|, are the consequence of oscillatory structure in the density of states of quarks and occur whenever 



the Landau level pass the quark Fermi surface 



14, 



2i. 



In Fig. mb), at T = 70 MeV and in the regime for eB < eBt, the diquark condensate is not built 
at all, i.e. A^ = 0. This result indicates that a normal phase exists at T > 70 MeV, n = 460 MeV and 
in the regime eB < 0.65 GeV'^, and that a second order phase transition from the normal to the CSC 
phase occurs at (T, ^) ~ (70,460) MeV and critical magnetic field eB ~ 0.65 GeV^. These results 
coincide with the observations in the complete T — eB phase diagram from Fig. [TWi) for fi = 460 
MeV. Indeed, our numerical computations show that for fi = 460 MeV, the critical temperatures in 
the whole range of eB £ [0, 0.6] GeV^ are smaller than 70 MeV. In eB = GeV^ and fi = 460 MeV, 
it is Tc — 68.5 MeV and decreases with eB < 0.5 GeV^ [see also Fig. HH^i)]. Using this value of Tc in 
zero magnetic field and the numerical value of the diquark gap at zero temperature and zero magnetic 
field from Fig. [5^b), i.e. Ab(T = 0) ~ 127.5 MeV, the above mentioned BCS ratio can be computed. 
It is given by 

Tc 68.5 

~ ~ 54 

Ab{T = 0) " 127.5 " ' 

which is in good agreement with the BCS ratio ^^^^q) = e"'ir~^ + 0{g) ~ 0.567 + 0{g) of critical 



temperature to the zero-temperature 2SC mass gap of QCD in zero magnetic field [37(]. Here, 7 ~ 0.577 
is the Euler-Mascheroni constant. As for the regime of strong magnetic fields, eB > 0.6 GeV^ at fixed 
T = 70 MeV, Ab increases monotonically with eB, as expected from piL12p for n = 460 MeV, and 
A as weh as Gd from piL2j) . 



See also the explanation in the paragraph following Fig. 1141 where we have compared Fig. [T4ji) with Fig. [6jb). 
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FIG. 7: (a) The dependence of the chiral gap ctb on eB for fixed T = 100 MeV, at various chemical potentials 
/X. (b) The dependence of diquark gap As on eB for fixed T = 50 MeV, at various chemical potentials. 

In Fig. [71 the dependence of chiral and diquark gaps on the external magnetic field for constant 
T are demonstrated for various fi. As it is shown in Fig. EH^a), the chiral mass gap aB decreases by 
increasing fi. Small van Alfven-de Haas oscillations occur in the regime below the threshold magnetic 
field eBt ~ 0.5 GeV'^ and die out in the linear regime above eBf. In this regime as monotonically rises 
with eB, and the effect of fj, is minimized. In the regime below eBt, for /u = 300 MeV, the chiral mass 
gap as vanishes. This indicates a first order phase transition from the normal into the xSB phase at 
fixed T = 100 MeV, fi = 300 MeV and for eB ~ 0.5 GeV^. This result agrees with our findings from 
the complete T — eB phase diagram plotted in Fig. [HTc) for fixed fj, = 300 MeV. As it can be checked 
in Fig. [Hlfc). at fixed T = 100 MeV, the xSB phase appears first for eB > 0.5 GeV^. What concerns 
the ei?-dependence of the diquark mass gap in Fig. [7|^b), it increases with fi in the regime below the 
threshold magnetic field eBf. In the regime above eBt, however, Ab decreases by increasing fi. Same 
phenomenon was also demonstrated in Fig. W^c) at low temperature T = 20 MeV and for eB = 0.5 
near eBt, and was shown to be in full agreement with the analytical result (1111.12^ [16]. Note that 
strong van Alfven-de Haas oscillations are responsible for vanishing in the regime eB G [0.4,0.6] 
GeV^. They induce CSC-Normal-CSC second order phase transitions, as it can be seen in the same 
regime of eB in Fig. HH^i). The plots in Fig. Ul^h) show also that the critical magnetic fields for the 
transition from the normal into the CSC phase increase by increasing the chemical potential [see Fig. 
[Ib)].i3 



Detailed analysis on the dependence of critical eB on /i and vice versa will be performed in Sec. III.B. 
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3. The T -dependence of a B and Ab 



In Fig. [HI the temperature dependence of the chiral condensate is presented for fixed /i = 250 MeV 
(panel a), fi = 300 MeV (panel b), and various magnetic fields eB = 0,0.3,0.5 GeV^. Whereas 
the transition from the xSB into the normal phase for /i = 250 MeV is of second order [continuous 
decreasing of as to as = in Fig. EJ^a)], it is of first order for ^ = 300 MeV [discontinuous transition 
from as ^ to as = in Fig. [8)^b)]. As it is also expected from our previous results, the magnetic 
field enhances the production of 0"^, that increases by increasing the value of the magnetic field. 



a)n = 250 MeV b) /j = 300 MeV 




T (MeV) T (MeV) 

FIG. 8: The T-dependence of <tb for fixed /i = 250 MeV (panel a) and /i — 300 MeV (panel b) for various 
eB = 0, 0.3, 0.5 GeV^. Whereas for fi — 250 MeV the transition from the chiral to the normal phase is of second 
order (blue curves in panel a), for fi = 300 MeV, this transition is of first order (green dots in panel b). 

The temperature dependence of is presented in Fig. [9]for fixed fj, = 480 MeV and eB = 0, 0.4, 0.5 
GeV^ (panel a) , below the threshold mag netic field eBt ~ 0.45-0.50 GeV^, as well as for fixed fi = 480 
MeV and eB = 0.6,0.7,0.8 GeV^, above the threshold magnetic field. As it turns out magnetic fields 
below the threshold magnetic field, suppress the production of Ab, whereas magnetic fields stronger 
than the threshold magnetic field enhances its production [see also Fig. [6jb)]. Similarly, for eB < eBt, 
the value of critical temperature decreases by increasing the magnetic field, whereas for eB > eBt, 
the critical temperatures increase by increasing the magnetic field. 

The same effect is also observed in Fig. [lOl In Fig. [TOl a). the T-dependence of ctb is plotted 
for constant eB = 0.3 GeV^ and various /i = 100, 200, 300 MeV. At low temperature T < 50 MeV, 
aB remains constant, and decreases by increasing the temperature. Increasing the chemical potential 
only accelerate this decrease, i.e. the critical temperature decreases by increasing /i. Whereas the 
transition at = 100, 200 MeV is of second order (blue dots), a first order phase transition occurs at 

The eB dependence of critical temperatures will be discussed in detail in III.B. 
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FIG. 9: The T-dependence of As for fixed ^ = 480 MeV for various eB = 0,0.3,0.5 GeV^ below tlie tlireshold 
magnetic field (panel a) and for eB — 0.6,0.7,0.8 GeV^ above the threshold magnetic field (panel b). Below 
(above) the threshold magnetic field the production of diquark is suppressed (enhanced) . 

/i = 300 MeV (red dots). The T-dependence of diquark condensate is demonstrated for eB = 0.3 
GeV^ below the threshold magnetic field [Fig. [ToT b)]. and for eB = 0.7 GeV^ above the threshold 
magnetic field [Fig. [TOlc)]. As it turns out, whereas the diquark condensate increases by increasing 
/_i for fixed T and eB < eBt [Fig. nUTb)]. it decreases by increasing /i for eB > eBt [Fig. [TUb]. This 
result coincides with our results from Fig. [7|^b). Moreover, as we can see in Figs. [ToTb) and 8(c), the 
diquark condensates decreases by increasing the temperature. Similar to what is observed in Fig. |9l 
for eB = 0.3 GeV^, below eBt, the critical temperature increases by increasing [Fig. [TOTb)]. while 
for eB = 0.7 GeV^ above eBt, the critical temperature decreases by increasing ^. Later, we will see 
that this result is also in full agreement to our result from Fig. [Hlfi). 



a)eB = 0.3 GeV^ b) eB = 0.3 GeV^ c) eB = 0.7 GeV^ 




T (MeV) T (MeV) T (MeV) 

FIG. 10: T-dependence of as and As for different for fixed eB. 
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B. Phase diagrams of hot magnetized 2SC quark matter 



To complete our study on the effect of (T, /j,, eB) on the quark matter including mesons and diquarks 
in the presence of constant magnetic field, we will present in this section the phase structure of the 
model in a T — /X plane for various fixed eB in Fig. [TTl The T — eB phase diagram of the model 
for various fixed is presented in Fig. [T^ and finally in Figs. [T6land [T7l the /i — eB phase diagram 
is explored for various fixed T. As in the previous section, green dashed lines denote the first order 
phase transitions and the blue solid lines the second order transitions. 

1. T — II phase diagram for various fixed eB 
In Fig. [TTl the T — fj, phase diagram of the model is demonstrated for various fixed eB. 



eB = OGeV^ eB = 0.1GeV^ eB = 0.2GeV' 




fi (MeV) II (MeV) /i (MeV) 



FIG. 11: The T — fi phase diagram of hot magnetized 2SC quark matter is presented for various eB. Blue 
solid lines denote the second order phase transitions and the green dashed lines the first order transitions. The 
critical and tricritical points are denoted by C and T, respectively. 

As we have argued at the beginning of this section, three different phases appear in this system; 
the xSB phase, the CSC phase and the normal phase. To study the effects of constant magnetic fields 
on the phase transition in the T — fi plane, let us compare the plots from Fig. [TT]with the plots from 
Figs. [2]m As we have observed in Figs. [2]^a)-[2]^c), at a fixed and low temperature T < Ti, first order 
transitions occur between the xSB and the CSC phase. By increasing the temperature to Ti < T < T2 
after a first order phase transition between the xSB to the normal phase, the normal phase goes over 
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TABLE I: Critical and tricritical points in T — /i phase diagrams of Fig. Ill) denoted by (/icrjTcr) and {fitr,Ttr), 
respectively. Here, eB is in GeV^, T and fi are in MeV. 

into the CSC phase in a second order phase transition [see 2(a)-2(b)]. At higher temperature T > T2, 
no CSC phase exists. Here, depending on the strength of the external magnetic fields, only second or 
first order transitions occur between the xSB and the normal phase [Figs. [4l^a)-IHc)]. As it turns out, 
the values of Ti and T2 depend on the external magnetic field, as it can be seen also in the plots from 
Fig. M 

The phase transition from xSB to CSC phase and from CSC to the normal phase are always of 
first and second order, respectively. The order of phase transition between the xSB and the normal 
phase depends, however, on the value of the external magnetic field eB. To describe this effect, let us 
consider the blue solid lines in Fig. [TTl that demonstrate the second order phase transitions between 
the xSB and the normal phases. Starting from high temperature and zero chemical potential, they all 
end at the critical points denoted by C (black bullets). They are then followed by first order critical 
lines (green dashed lines) between the xSB and the normal phase. The latter start at the critical 
points C and end up at the tricritical points T (red bullets), where three phases coexist. In Table [H 
the values of the temperature and chemical potentials corresponding to the critical points (Tcr.,/Xcr-) 
and the tricritical points (Ttr,fJ^tr) are presented. 

As we can see in Fig. [TTl by increasing the external magnetic field, the black bullets are shifted 
more and more to higher values of temperature and chemical potential (see also second column in 
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Table H]), so that the distance between the black and red bullets increases by increasing the strength of 
external magnetic field. This is another effect of the external magnetic field; increasing the magnetic 
field strength changes the type of the xSB to normal phase transition from second to first order. This 
can also be observed in Figs. [ll^a)-|3I^c), where at high enough temperature, the xSB phase transition 
to the normal phase is first of second order and then, for larger magnetic fields, changes to a first 
order phase transition. 

Increasing the magnetic field strength leads also to an increase in the values of {^tr^Ttr) corre- 
sponding to the tricritical point (see the third column of Table H]). This confirms the conclusion at the 
end of Sec. III. A. 3, where it was stated that the magnetic field above a certain threshold magnetic 
field enhances the production of the diquark condensate A^, and the CSC phase can therefore exists 
up to T ~ 100 MeV. In what follows, we will show that above certain threshold magnetic field, eBt, 
the analytical results arising in a LLL approximation are comparable with the numerical results in- 
cluding the contributions of all Landau levels. To do this, we will compare the analytical expression 
for the second order critical lines of a transition between the xSB and the normal phase with the 
corresponding numerical data to determine the threshold magnetic field for the LLL approximation. 
Similar comparison will be then performed for the second order critical line of the transition between 
the CSC and the normal phase. 



eB = 0.5 GeV^ eB = 0.7 GeV^ 
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FIG. 12: The second order critical lines of transition between the xSB and the normal phase. The red dots 
denote the analytical data that arise in a LLL approximation and the blue solid line the numerical data including 
the contributions of all Landau levels. For eB = 0.7 GeV^ stronger than eBt, the analytical and numerical data 
exactly coincide. 



The second order critical surface of the transition between the xSB and the normal phase is de- 
termined explicitly in App. lA 1[ Here, we will present only the final results. Following the method 
presented in [2^, in the phase space of the thermodynamical parameters (T, /i, ei?), the second order 
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critical surface is determined by solving 
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2l| 



9aff(a,A = 0) 
lim —5 = 0. 

0-2^0 oa'^ 



(III.13) 



Setting n = in from (|II.23p and plugging the resulting expression in (IIII.13P , leads to the second 
order critical surface in the phase space of these parameters [see (jA.lOp in App. lA Ij 



1 



AGs 47r2 



dz [2z + j F[z; T, fi] = 0, 



with 



F[z-T,fj] 



sinh(/3z) 



cosh(/32;) + cosh(/3/i) 
iers of K = ^ and keej 
we arrive after a straightforward computation at [see also ()A.17p in App. lA Ij 



(III. 14) 
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Expanding the integral in piI.14p in the orders of k = ^ and keeping only terms of the order 
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+ 3eB 



dz 



tanh z 



where A = and a is defined by 



a{T, eB; A) = (tanh A + A tanh^ A - A) + 



3eB tanh^ A 3eB 
+ 



8r2 A 8T^ 
Moreover, Li2(2;) in piI.16P is the dilogarithm function defined by 

ln(l - z) 



dz 



tanh^ z 



U2{z) = 



dz 



(III. 16) 



(III.17) 



(III.18) 



To determine the second order critical line between the xSB and the normal phase in the T — fj, plane, 
we have to fix eB. The analytical results arising from the LLL approximation are demonstrated in Fig. 
[T2] by red dots. Blue solid lines denote numerical results for second order critical lines including the 
contributions of all Landau levels. In Fig. \T2\ the analytical and the numerical results for eB = 0.5 
GeV^, approximately at the threshold, and for eB = 0.7 GeV^, above the threshold magnetic field, 
are compared. The qualitative behavior of two curves coincides above the threshold magnetic field 
eBt — 0.5 GeV^, where the system is in the regime of LLL dominance. 

Similarly, the second order critical surface between the CSC and the normal phase is determined 



from 



20I, 



2ii 



9aff(^ = 0,A) 
hm = 0. 

A2^o dA^ 



(III. 19) 



Setting n = in ilea from pi.23P and plugging the resulting expression in piI.lQp . we arrive at [see 
(|A:2T]1 in App. [A2] 



eB~\T,f^;A) = ^ 



H 



A+ji 
2T 



+ H 



A 



2T 



(III.20) 
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where H{z) is 



HIz) 



E 

n=l 



(2n - l)(2n) 



--Br, 



(111.21) 



Here, i?„ are the BernoulH's numbers. Equation pil.20p leads to a relation between the phase space 
parameters {T,iJ,,eB). In Fig. [T31 we have fixed eB to be eB = 0.5,0.7 GeV^, and compared the 
analytical data of the second order critical line (red dots), arising from pil.20p - pil.2ip with the 
numerical data (blue solid line). For eB > eBt ~ 0.5 GeV^, the analytical and numerical data exactly 
coincide (see the plot of eB = 0.7 GeV^ in Fig. [T3l) . We conclude therefore that for eB > eBt the 
higher Landau level are decoupled and the properties of the phase transitions are essentially affected 
by the contribution from the LLL. 
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FIG. 13: The second order critical lines of transition between the CSC and the normal phase. The red dots 
denote the analytical data that arise in a LLL approximation and the Blue solid line the numerical data including 
the contributions of all Landau levels. For eB = 0.7 GeV^ stronger than eBt, the analytical and numerical data 
exactly coincide. 



2. T — eB phase diagram for various fixed /x 

In Fig. the T — eB phase diagram of hot 2SC quark matter is presented for various /u. Blue solid 
lines denote the second order phase transitions and the green dashed lines the first order transitions. 
The critical points are denoted by C and the tricritical points by T. The exact numerical values of 
the critical and tricritical points are presented in Table HIl 

In Fig. [IHa) andMlh), the critical lines are plotted for ^ = 0, 200, 250 MeV and = 270, 280 
MeV, respectively. For relatively small values of /_i < 270 MeV the transition between the xSB and 
the normal phase is of second order for the whole range of magnetic field eB £ [0,0.8] GeV^. When 
we increase the chemical potential to /i = 270, 280 MeV in Fig. HD^b), the type of the phase transition 
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FIG. 14: The T — eB phase diagram of hot magnetized 2SC quark matter is presented for various ^. Blue 
solid lines denote the second order phase transitions and the green dashed lines the first order transitions. The 
critical points are denoted by C and the tricritical points by T. 



changes from second to first order in the regime eB G [0.1,0.5] GeV^, below the threshold magnetic 
field. By increasing the magnetic field and entering the regime of LLL dominance, the first order phase 
transition goes over into a second order transition [see in Table |TI] for the exact values of the critical 

n n 

points]. The same effect has been previously observed in [21[| (see Fig. 4(b) in [21[ and compare it 



with Fig. [T^a) and I14lfb) of the present paper). In Fig. [Hlfc) the chemical potential is increased to 
fi = 300 MeV. Here, a small region of CSC phase appears in the regime eB G [0.1,0.3] GeV^, between 
two regions of the xSB phase. The critical points appearing in Fig. [HTa) are shifted to higher (T, eB) 
[see also Table Hi] . The small CSC "island" enlarges by increasing the chemical potential to ^ = 310 
MeV in Fig. [Hlfd). and remove the xSB phase appearing at eB < 0.15 GeV^ totally at ^ = 325 MeV 
[Fig. fl^e)]. This is expected from Figs. [2]J^a)-l2]J^c), where the critical chemical potential from the 
transition of xSB to CSC phase is approximately ~ 325 MeV. The xSB phase appearing only for 
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T — eB phase diagram (Fig. I14p 

^ {eBcr,Tcr)l {eBcr,Tcr)2 {eBtr,Ttr)l {eBtr,Ttr)2 

< 250 - 

270 (0.16,97) (0.310,128) 

280 (0.10,83) (0.510,159) 

300 - (0.578,190) (0.165,32) (0.290,39) 

310 - (0.595,190) (0.095,43) (0.310,44) 

325 - (0.635,209) - (0.485,50) 

340 - (0.668,221) - (0.520,54) 

370 - (0.740,250) - (0.585,64) 

> 400 - 



TABLE II; Critical and tricritical points in eB — T phase diagrams of Fig. [TH denoted by {eBcr,Tcr)i,2 and 
{eBtr,Ttr) 1,2, respectively. Here, eB is in GeV^, T and fi are in MeV. 

strong magnetic field survives, but it is puslied away to the regime of LLL dominance by increasing 
the chemical potential to fi = 340 and /u = 370 MeV in Figs. fTUf f) and fTllf g). We notice that the 
transition from xSB to the CSC phase is always of first order, whereas the type of the phase transition 
from xSB to the normal phase depends on the external magnetic field. This confirm our conclusion 
from the previous section (See Fig. [IT] and the discussion in Sec. III.B.l). Increasing the chemical 
potential to /u > 400 MeV, the xSB region is completely removed from the interval eB £ [0, 0.8] GeV^, 
and a second order phase transition appears between the CSC and the normal phase in Fig. [Hlfh). In 
Fig. Hm^i), the second order phase transition is plotted for /x = 460, 500 MeV.^^ Above the threshold 
magnetic field, the critical temperature increases by increasing the magnetic field. This may open the 
possibility to observe the 2SC phase in future heavy ion experiments. 

Let us also emphasize that the plots in Fig. [THare in very good agreement with plots from previous 
sections. As an example, let us consider Fig. [6jb), the eB dependence of at fixed /i = 460 MeV, 
and compare it with the curve // = 460 MeV in Fig. [T4lfi). Let us then focus on the range of 
eB G [0.4,0.6] GeV^ in both figures. As it turns out from Fig. EJ^b), at low temperature T = 20 
MeV, Ab ~ 70 - 80 MeV in the interval eB G [0.4, 0.6] GeV^, while at higher temperature T = 60, 70 
MeV, Ab vanishes and apparently a normal phase appears in this regime of eB. The appearance of a 
normal phase can be checked in Fig. [T^i). For fixed T = 20 MeV, there is a CSC phase in the regime 
eB G [0.4,0.6] GeV^, while at higher temperature T = 60 MeV and T = 70 MeV, a normal phase 

^'^ In Fig. HUi), we have changed the scale of the plot from T e [0, 360] MeV to T £ [35, 125] MeV, in order to magnify 
the van Alfven-de Haas oscillations that appear in the regime below the threshold magnetic field eBt — 0.5 GeV'^. 
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appears in the same regime of eB. This phenomenon is because of van Alfven-de Haas oscillations 
that occur in this regime of eB £ [0.4, 0.6] GeV^ in Fig. [6] as well as in the phase diagram of Fig. 
[T^i).^^ Moreover, the fact that increases for eB stronger than a certain threshold magnetic field 
eBf shows the important interplay between the effect of temperature and external magnetic field on 
the formation of bound states, which could be potentially relevant in connection with the question of 
accessibility of 2SC phase in experiments. Whereas below the threshold magnetic field, by increasing 
the temperature, the diquark condensate is destroyed and a normal phase is build up, the production 
of As is enhanced by magnetic fields stronger than eBt in a corresponding CSC phase. According to 



our previous studies in [16(], we believe that above eBt the system enters the LLL dominant regime, 
where the effect of magnetic catalysis enhances the formation of mass gaps.^^ 

To show how the numerical results including the contributions of all Landau levels coincide with the 
analytical results consisting only of the contribution of the LLL above the above mentioned threshold 
magnetic field, we consider, as in the previous section, equation piL14p . expressing the second order 
critical surface between the xSB and the normal phase. Setting for instance ;U = 0, we get, as it is 
shown in App. lA 11 the second order critical line in the T — eB phase [see ()A.12p ] 

4ir^ f 1 A' T'' 



Li2(-e-2^) -2Aln(l + e-2^) 



(IIL22) 



where A = A/2T. Moreover, the dilogarithm function Li2(2;) and ii{z) are defined in piLlSp and 
piL2ip . respectively. In Fig. [T5l the analytical and numerical data are compared. Red dots denote 
the analytical data and the blue solid line the numerical data. Whereas at the threshold magnetic 
field eBt ~ 0.5 GeV^ the qualitative behavior of both data are similar, for eB > eBt the analytical 
and numerical data exactly coincides. 



3. /i ~ eB phase diagram for various fixed T 

In Figs. [16] and [T7\ the ^ — eB phase diagram of hot magnetized 2SC quark matter is illustrated 

for various fixed temperature T = 20, • • • 200 MeV. Green dashed lines denote the first order phase 

transitions and blue solid lines the second order phase transitions. The critical points are denoted by 

C (black bullets) and the tricritical points by T (red bullets). 

Let us consider the first plot in Fig. [161 the u — eB phase diagram at T = 20 MeV. This plot 

I I 

is similar to the phase diagram, which was found in ^] at T = 0. It includes a first order phase 



Similar van Alfven-de Haas oscillations occur in the "gap vs. eB / graphs" of a magnetized CFL model studied in 
[l3 | and [13] • In [l^, where the effect of magnetic fields on the 2SC gap is studied at T = 0, same oscillations occur 
in the same regime of eB for /x = 460 MeV (see Fig. |3K in ,16|). 

See the explanation in Sec. I, for different mechanisms being responsible for the production enhancement of meson 
and diquark mass gaps by strong magnetic fields. 
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FIG. 15: The second order critical lines of transition between the xSB and the normal phase. Red dots denote 
the analytical data that arise in a LLL approximation and the blue solid line denotes the numerical data including 
the contributions of all Landau levels. The threshold magnetic field is at eB ~ 0.5 GeV^. For eB > eBt the 
analytical and numerical data exactly coincides. 

transition between the xSB and the CSC phase. The latter goes over into the normal phase in a second 
order phase transition. This plot confirms our findings in Figs. [2]^a)-[2]^c) and Fig. [TTJ At T ~ 50 
MeV, the CSC phase is broken into two separated islands by the normal phase, and a tricritical point 
(red bullet) appears at {eB, /i) = (0.490 GeV^, 333 MeV). The normal phase between two CSC islands 
appears at eB £ [0.4,0.6] GeV^. This confirms our results from[71^b), where in the same regime of eB, 
the CSC condensate vanishes and a normal phase occurs. By increasing the temperature both 
CSC islands shrink, so that at T = 70 MeV only three separated CSC islands are remained in the 
regime below the threshold magnetic field eB < 0.5 GeV^. They are then totally destroyed at T > 100 
MeV and in the relevant interval eB £ [0,0.8] GeV^. Moreover, by increasing the temperature from 
r ~ 50 MeV to T ~ 100 MeV, the CSC island and the corresponding tricritical point that exist in 
the regime above the threshold magnetic field are shifted away to higher values of fi and eB (for the 
exact values of the tricritical points see Table . 

At T > 100 MeV, the fj, — eB phase space consists only of the xSB phase and the normal phase. 
The critical points that appear at T = 100 MeV at (0.14 GeV^ 267.3 MeV) is shifted to higher values 
of eB and fi as it can be seen in Fig. [T7] and Table IIIII High temperature suppresses the production 
of the chiral condensate in the regime of small magnetic fields (compare the curves in Fig. [17] in the 
regime eB £ [0,0.4] GeV^). This effect is compensated by increasing the magnetic field to above the 

As we have explained in Footnote 2, in a three-flavor NJL model at low temperature and fj, > 500 MeV, the CSC 
phase goes over into a CFL color-superconducting phase. In our model, where no CFL phase can be built, the second 
order transition between the CSC and the normal phase at low temperature and fi > 500 MeV is only assumed to 
exist. The same assumption is also done in [l^ . 
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FIG. 16: The /i — eB phase diagram of hot magnetized 2SC quark matter is presented for various T. Blue sohd 
hues denote second order phase transitions and the green dashed hues first order transitions. The critical points 
are denoted by C and the tricritical points by T. The brown regions denote the CSC phases. A small CSC 
phase appears at T = 100 MeV on the right side of the last plot. The tricritical point (red bullet) is shifted to 
the regime eB > 0.8 GeV^ (see Table Hill . Black bullet denotes the critical point. 

threshold magnetic field. 




FIG. 17: fi-eB phase diagram for T = 120, 150, 180, 200 MeV (from top to bottom). 
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fi — eB phase diagrams (Figs. I16L I17p 



T 




{eBtr, fJ.tr) 


^ 20 


— 
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50 


— 


(0.490,333) 


60 


— 


(0.554,356) 


65 


— 


(0.580,365) 


70 


— 


(0.610,387) 


100 


(0.140,267) 


(0.774,467) 


120 


(0.240,245) 




150 


(0.458,269) 




180 


(0.520,273) 




200 


(0.610,314) 





TABLE III: Critical and tricritical points in eB — /i phase diagrams of Figs. \W\ and [T71 denoted by {eBcrjfJcr) 
and {eBtr, fJ.tr), respectively. Here, eB is in GeV^, T and are in MeV. 

IV. CONCLUDING REMARKS 



In this paper, we have introduced a two-flavor superconducting NJL type model including mesons and 
diquarks at finite temperature T, chemical potential fi, and in the presence of a constant (rotated) 
magnetic field eB. We were in particular interested on the effect of (T, fi, eB) on the formation of 
chiral and color symmetry breaking bound states, and on the nature of phase transitions. One of the 
most important effects of the magnetic field is its competition with T and fi in the formation of bound 
states. This effect is demonstrated in Figs. [6] and [71 As it is illustrated in Fig. [6l although it is 
expected that increasing T suppresses the formation of chiral and color condensates, it turns out that 
this specific effect of T is minimized in the regime of strong magnetic field, for eB larger than a certain 
threshold magnetic field eBt — 0.5 GeV^. However, whereas the chiral condensate as increases by 
increasing eB in the whole range of eB E [0,0.8] GeV^ [Fig- Efa)]) there are regions below eBt, where 
the color condensate vanishes [Fig. E^b)]. Same phenomenon is also observed in Fig. \7\ where, 
at fixed temperature and various chemical potentials, the dependence of aB and on the magnetic 
field is plotted. Here, as in the previous case, the effect of fi in suppressing the production of the chiral 
and (2SC) color condensates is compensated by eB in the regime above the threshold magnetic field. 

According to our results from [ig], we believe that in the regime above eBt, the dynamics of the 
system is solely dominated by the lowest Landau level. In 16|, a magnetized 2SC model at zero 
temperature is considered, and, the numerical value of the threshold magnetic field is determined 
by comparing the analytical eB dependence of cr^ and A^, arising in a LLL approximation, with 
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the corresponding numerical data including the contribution of all Landau levels. According to this 
comparison, the threshold magnetic field at T = turned out to be eBt ~ 0.5 GeV^. In the present 
paper, the same model is considered at finite temperature. Comparing the analytical expression of 
the second order critical line of the transitions from xSB to the normal phase, arising from a LLL 
approximation, with the corresponding numerical data including the contribution of all Landau levels 
in Fig. [151 we arrive at the same threshold magnetic field as in T = case. Note that in {igI . the 
color neutrality condition was imposed on the magnetized 2SC quark matter. A comparison between 
the p.esented .suits .n th.s pap., where no color neut»ht, ts ..posed, and the results frou, Q. 
shows that color neutrality, being very small [35l|, has no significant effect on the value of the above 
mentioned threshold magnetic field. 

The nature of the phase transition is also affected by the external magnetic field. The most 
significant effect in this regard is illustrated in Fig. [TTl where the T — jj, phase space is plotted for 
various eB. As it is demonstrated in this figure, the phase space of the system includes three different 
phases. We observe that the distance between the black and red bullets, denoting the critical and 
the tricritical points, increases by increasing the strength of the external magnetic field. This implies 
that by increasing eB and keeping (T, ^u) fixed, the second order transition from the xSB phase to 
the normal quark matter changes into a first order transition. This can also be observed e.g. in Fig 
[TTl on each isothermal critical curve, a second order phase transition changes into a first order one by 
increasing the strength of the magnetic field. Similar effect occurs in the electroweak phase transition 
by implying an external hypermagnetic field ^ and in type I QED superconductivity. Whereas the 
type of the phase transition from xSB to the normal phase changes by increasing the strength of eB, 
the latter has no effect on the nature of the phase transition between the xSB and the CSC phase 
(first order) and between the CSC and the normal phase (second order) [see the phase diagrams in 
Figs. [TT| [T4l I16j . The only crucial effect of eB concerning the CSC phase, are the observed CSC- 
Normal-CSC second order phase transition induced by strong van Alfven-de Haas oscillations in the 
regime eB G [0.4,0.6] GeV^ [see e.g. Figs. EJ^b) and 1141^1) at fixed T = 60 MeV], and an increase of 
the tricritical temperature and chemical potential by increasing the magnetic field. The latter can be 
observed in Table |T] for eB > eBt ~ 0.5 GeV^. This could be relevant in relation with the question 
addressin g th e accessibility of 2SC phase in present or future heavy ion experiments, which is recently 
32] . We have extended the setup of the model used in 



posed in 



32] by considering the effect of 



constant magnetic fields, which are supposed to be created in the non-central heavy ion collisions 
and are estimated to be of order eB ~ 0.03 GeV^ at RHIC energies, and eB ~ 0.3 GeV^ at LHC 



See Footnote 9. 
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energies. These amounts of magnetic fields are, according to our observation in this extended model, 
far below the range of magnetic fields which could have significant effects on the CSC phase transition 
by partly compensating the effects of {T,fi), as described above. However, we believe, that more 
realistic models are to be examined to find a satisfactory answer to this interesting question. 

The model which is used in this paper can be extended in many ways, e.g, by considering the 
2SC-CFL phase including color neutrality. Similar computation can also be performed within a color 
neutral PNJL model. It is also intriguing to explore the effect of the axial anomaly in the same 
context. 
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Appendix A: Second order critical lines in the LLL approximation 

In this appendix, we will determine the main equations leading to the second order critical lines cor- 
responding to the transition between the xSB and the normal phase, as well as the transition between 
the xSB and the CSC phase, in the limit of very strong magnetic field in the LLL approximation. 



1. Transition between the x^B and the normal phase 



(T,u,eB), the second order 



In the phase space spanned by the intensive thermodynamical parameters 
critical surface between the xSB and the normal phase is determined by 

I™ 7^ = 0- 

Using Q^ff from ()II.23p and setting n = 0, to consider only the contribution of the lowest Landau level, 
we get first 



(A.l) 



cP'p f°° ds p + Kfi 



AGs 2V^^^^\J (27r)3 7o ^ p 

x[l + 2^(-l/ e-T^] =0. 
t=i 

Denoting the temperature independent part of ()A.2p with It=Oi and using 

1 ds X- K/i ^_,(^_^^)2 ^ 9{x- n) 



+seB r — r '^^^ P3+^/i ^-.few 

Jo Vs Jo 47r2 p3 

(A.2) 



(A.3) 
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from [38|, we get 



It=o = -TT^- 7^ 9{p -fi)- 3eB 



4Gs Jo 27r2 7o 47r2 

1 (A2-/i2) 3eS^ /A 



In - . (A.4) 



AGs 47r2 47r2 \fJ.J ' 

To evaluate the momentum integration, we have introduced the ultraviolet cutoff A. Note that (|A.4p 
is comparable with our results from [la]. The temperature dependent part of (|A.2p . which is denoted 
by It^O) can be evaluated using 

^ jT" = VI (e-/3^l-Mlsgn(z - ^) + e-^'(^+^^) , (A.5) 

and is given by 

It^o = f;(-l/{^'|^e-^'^'^sinh(^/3z)(^ + 2z)-|^|^e-^^^cosh(W (^ + 2z)}, 

(A.6) 

where z is a generic integration variable replacing p = |p| or p3 in the one dimensional integrations of 
()A.2p . To perform the summation over i, we use 

^(-l/sinh(^/3z)e-^'^'^ = \F[z;T,fi], 
1=1 

oo 

^^(-l)^ cosh {IM e-'P^ = - {F[z; T, /.]-!), (A.7) 
£=1 

where 

1 smh.(Bz) 
F[z;T,f,]^ ^ (A.8) 

cosh(pzj + cosh(p/ij 

Plugging (|A.7p in ()A.6P and performing the integration over z as far as possible, the temperature 
dependent part It^o is given by 

The second order critical surface of the transition between the xSB and the normal phase as a function 
of the phase space variables (T, /x, eB) is then given by adding (jA.4p and (|A.9p and reads 

1 1 f 3gB \ 

-—^ dz[2z + F[z;r,^]=0, (A.IO) 



AGs Att"^ Jq \ z J' 

where F[z;T, fi] is defined in (jA.Sp . Note that ()A.10p is only valid in the limit of strong magnetic field, 
where a LLL approximation is justified. Choosing /i = 0, (jA.lOp reads 

AGs 



1 1 /"^ , 3eS\ Bz 
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This leads to the fohowing relation between T and eB 



47r2 J 1 



A2 



Li2(-e~^^) -2Aln(l + e 



-2X\ 



(A.12) 



where A = A/2T and 



H{z) = ^ .~ ^ 



n=l 



(2n - l)(2n)! 



(A.13) 



and Li2(z) is the dilogarithm function defined by 



Li2(2;) = - I dy 



ln(l - y) 



(A.14) 



/o y 

and Bn are the Bernoulli's numbers. This result will be used in Sec. III.B to determine the threshold 
magnetic field for the LLL approximation at /i = 0. For arbitrary values of the chemical potential, 
we expand (|A.10p in the orders of ^ and keep terms up to C'((^)^), and arrive at the explicit (T, //) 
dependence of eB 



eB{T,fi;A) 



1 



IT 



37 I Gs 



A^-Li2(-e-2A)^2Aln(l + e 



0-2A 



12 



-11^ (tanli A + A tanh^ A — A) | , 



where 7 is defined by 



fi"^ tanh A 



8T2 A JO 
Similarly, the (T, eB) dependence of /x is given by 



^ ( tanh z tanh^ z 

dz 



8r2 z2 



(A.15) 



(A.16) 



^\T,eB;A)^-\-^ + iT^ 
a \ Gs 



A2-Li2(-e-2^) + 2Aln (l + e" 



2A^ _ !L 
12 



A 



tanh z 



+ 3eB / dz- 

'0 z 



(A.17) 



where a is defined by 



(T, e5; A) = (tanh A + A tanh^ A - A) + 



3eB tanh^ A 3eB 
+ 



dz- 



tanh^ z 



(A.18) 



8r2 A 8T2 7o 

Relations (|A.17p will be evaluated numerically in Sec. III.B to determine the threshold magnetic field 
for the LLL approximation. 



2. Transition between the CSC and the normal phase 



For the transition between the CSC and the normal phase, we use the equation 20|, |21[ 



lim TT-Tn = 0. 

A2^o aA2 



(A.19) 
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Setting n = in from pi.23p and plugging the resulting expression in ()A.19p . we get 



eB 



K = ±l 



HID j ds j aps 



dP3 



Jo 47r2 



(A.20) 



Evaluating the T-dependent and independent parts of the expression in the r.h.s. of (IA.20p using (lA.Sp 
and (jA.Sp . respectively, we arrive after some straightforward manipulations at the relation between 
the phase space parameters (T, /i, eS), 



G 



TT' 

Gd 



2T tanhz 
dz h 



z 



~2T- tanhz 
dz 



z 



H 



A + /X 
2T 



+ H 



A - /u 
2T 



(A.21) 



where H{z) is defined in (jA.lSp . Fixing one of these parameters in ()A.2ip . the second order critical lines 
for the transition from CSC to the normal phase arises in the phase space of two other parameters. 
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